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Abstract 

We study nuclear shadowing in J/ip and T production in hadron-nucleus 
interactions and in nucleus-nucleus collisions at the Relativistic Heavy Ion 



^ ■ Collider and the Large Hadron Collider. As a consequence of the perturba- 

tive Q 2 -dependence of gluon shadowing, we predict that T production is less 
suppressed than the J/ip. We show that antishadowing leads to enhanced 
J ftp production at xj < 0, an effect reduced for T production. 
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The production of charmonium and bottomonium states in hadron-nucleus interactions 
does not increase linearly with the nuclear target size A. Additionally, the ratio of hadron- 
nucleus to hadron-nucleon production, Sa = VhA/Auhp, is not constant but decreases as 
a function of the fraction of the center of mass momentum, Xf, carried by the produced 
resonance. There are several effects that contribute both to the A dependence of the total 
cross section and to the ratio Sa including nuclear absorption [|IJ, comover interactions |§, 
intrinsic heavy quark states ||, projectile energy loss |§, and nuclear shadowing fUJ. 

In this paper, we examine the role of shadowing of the nuclear parton distributions in 
quarkonium production. Nuclear shadowing modifies the target parton distributions so that 
xq A (x) 7^ Axq N (x). Deep inelastic scattering (DIS) data on the ratio Ra = F^/F® |7||| 
show that: i) nuclear shadowing begins to set in at xbj ~ 0.06 for all nuclei (Ra < 1); ii) 
Ra has a very weak Q 2 dependence; and Hi) an enhancement, or antishadowing, exists for 
0.06 ;$ xbj ~ 0.25 (Ra > 1). In addition, the NMC data on J/if) production with Sn and 
C targets show an enhancement of (13 ± 8)% for 0.05 < xbj ~ 0.3 ||, interpreted as gluon 
antishadowing. 

The weak Q 2 dependence of Ra is consistent with the models of JTU| , |rT| |, where nuclear 
shadowing is considered as a leading-order QCD effect with perturbatively generated Q 2 
dependence. These models differ in the treatment of the nonperturbative contribution at 
the initial scale Q^. We adopt the QCD Aligned- Jet Model (QAJM) |T| which accounts for 
the DIS data [p) , |P2]1 to study the effect of nuclear shadowing on Sa- 

We would like to point out two major consequences of shadowing and antishadowing 
for J/ip and T hadropro duct ion based on very general arguments. Since mj is three times 
larger than m^, the T is primarily produced within the antishadowing region of the target 
at present energies. However, the parton distributions are evaluated at Q 2 = m\, thus Q 2 
evolution reduces the effect of both shadowing and antishadowing in T production with 
respect to J/ip. The antishadowing apparently present in both the DIS measurements of Ra 
|7|,|3|] and the gluon distribution || reflects on Sa- We include the observed antishadowing 
effect, a general consequence of baryon number and momentum conservation in a nucleus, 



and explore the consequences of Q 2 evolution on Sa- 

We do not expect our shadowing model to fully explain the shape and magnitude of 
Sa- To quantitatively show the effects of nuclear shadowing, we have not included other 
A dependent effects. Nuclear absorption and comover interactions do not have a strong Xf 
dependence and mainly affect the magnitude of Sa [j],@] • O n the other hand, intrinsic heavy 
quarks |§ and projectile energy loss f§] have a significant Xf dependence. 

The leading order QQ bound state cross section is the integral of the free QQ production 

cross section from the QQ threshold to the meson pair threshold, 

dor p2m H /^ H vt (xi,x 2 ;x 1 x 2 s) , , 

- IF \ rdr pK il- 2 ' ; , (1) 



dx f JTmQ/y/i ^ x ) + 4r 2 

where m c = 1.5 GeV, m h = 4.75 GeV, m D = 1.867 GeV, and mg = 5.28 GeV. The bound 
state fraction, F, is a parameter that cancels in Sa- We include quark-antiquark annihilation 
and gluon-gluon fusion subprocesses in the convolution formula H pt , 



,7=> ™2> 



H pt (xi,X2-, m ) = G p (x 1 )G t (x 2 )(r(gg -> QQ; m 



J2 (l P (xi)q t (x2) + q P (xi)q t (x 2 ))cx(qq -+ QQ; m 2 ) , (2) 

q=u,d,s 



where x% t 2 = \{±Xf + Jx 2 + 4r 2 ) are the projectile and target momentum fractions. We 



use the MRS D— ' |13| distributions that describe the recent HERA data on F 2 P [14 . 

Next-to-leading order contributions to the free QQ cross section have been calculated. 
The Xf distribution is approximately that from lowest order multiplied by a theoretical K 
factor, cr(ttg) /cr(a 2 ) |TJ|], assumed to be absorbed in F. At lowest order, the J/ip is produced 
with pt = since the cc pair is back-to-back. A p T dependence may be introduced through 
an intrinsic transverse momenta, qx, of the initial partons as we show later. For p? 3> m ves , 



higher-order corrections are needed ||16| . Since we concentrate on the Xf distributions and 
the A dependence, where pt ~ rn res is dominant, our calculation is lowest order. 

When the target is a nucleus, the parton distributions, G t (x 2 ), qt(x 2 ) and qt{x 2 ) in Eq. 
(0) are the nuclear medium modified distributions evaluated at Q 2 ~ vnL and Q 2 ~ m\ for 



J /if) and T production. We calculate the shadowed distributions at the initial scale, Q 
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Ql ~ 5 GeV 2 , and x < x s h ~ 0.04 according to the QAJM ]TT[ and evolve in Q 2 using the 
Dokshitser-Gribov-Lipatov-Altarelli-Parisi equations (DGLAP) [|Tj]]. Recombination effects 
within a single nucleon are negligible at x ^ 10~ 3 [|5|,|T8||. In the QAJM, shadowing is 
generated by the interaction of hadronic configurations originating from fluctuations of the 
virtual photon with small transverse momentum, or jets aligned along the direction of the 
virtual photon, q/|q|, with k± < ko± « 0.4 GeV. The transverse separation of the aligned 
jets is correspondingly large (~ 1 fm) and their phase space is restricted by a factor oc 
k\^jM 2 . In QCD the transverse momentum of the fluctuation can be large so that nonaligned 
hadronic configurations are simultaneously present. Such configurations interact differently 
in the nuclear medium. In Refs. ]ll],|l3|] h was assumed that aligned jets with large transverse 
separation interact in the nucleus as vector mesons. Nonaligned jets, with k± > ko± and 
a small transverse separation have a correspondingly small interaction cross section, giving 
rise to color transparency. 

The nuclear baryon number and momentum sum rules: 

\ [ dxV A {x, Q 2 ) = [ l dxV N (x, Q 2 ) , (3a) 

A Jo Jo 

\ [ dxx[V A (x, Q 2 ) + G A (x, Q 2 ) + S A (x, Q 2 )} = 
A Jo 

f 1 dxx[V N (x, Q 2 ) + G N (x, Q 2 ) + S N (x, Q 2 )} , (3b) 

JO 

must be fulfilled at every Q 2 . In Eq. (|3b| ) Vn(A) = uj + ^t ', S^(A) = 2(«t + d t + s t ), and 
Gjv(a) — G t , are the nucleon, N, and nuclear, A, valence, sea and gluon distributions. In 
order to satisfy the sum rules and simultaneously allow for shadowing at x < 0.04, the parton 
distributions must be enhanced at higher values of x. We assume that the enhancement is 
shared by the valence quarks and gluons since the sea quark enhancement is consistent 
with zero |L9| (see also Ref. |T2[). Recent data from E665 || show that the enhancement 
in R A is concentrated in the interval 0.06 ^$ x ;$ 0.25 and is A-independent. These two 
crossover points, assuming no sea enhancement, are sufficient to constrain the valence and 
gluon antishadowing. 

Gluon shadowing should be larger than sea quark shadowing since the interaction of 



a qqg configuration in the 7* wavefunction with the target should be stronger than a qq 
configuration [[2I],[nj]. Indeed, the perturbative interaction cross section of a small color 
octet with a nucleon is 9/4 larger than a color triplet of the same transverse size. Assuming 
a smooth connection between the perturbative and nonperturbative domains, <7 7 *_» g g fli jv = 
|cr 7 *_ >g g ) jv, leading to larger effects on the gluon relative to the sea and valence quarks [|20 



The overall enhancement found in |1^J20|] is consistent with the data of ||, used to constrain 
the shape of the gluon distribution in the antishadowing region. A weak Q 2 dependence for 
Ra is obtained, in accord with the data of @||, while a stronger Q 2 dependence is found 
for the gluon ratio, Rj[ = Ga/Gn- Similar results were also obtained in [^]. 

In Fig. 1 we show x 2 as a function of Xf for J/ip (solid curve), T (dashed) production for 
several energies to point out the shadowing and antishadowing regions. At 200 GeV, 1(a), 
the J/ip is antishadowed for ^ Xf ^ 0.6. At 800 GeV, 1(b), where data exists for J/ip and 



T from —0.2 < Xf < 0.65 f23fl , the J/ip lies in the shadowing region for Xf > but enters the 
antishadowing region at backward values. The T is mainly produced in the antishadowing 
region, but when Xf < 0, it is in the EMC region. Since the parton distributions are evaluated 
at m\ and the gluon distribution evolves faster than the quarks, we expect S\ > S\ for 
Xf > 0. In 1(c) and 1(d), the Xf range corresponds to —0.35 < y < 2.5 and \y\ < 1, the 
rapidity coverage of the PHENIX fH and ALICE p| detectors at RHIC (y/s = 200 A GeV) 



and the LHC (\/i = 5.5 TeV), respectively. At RHIC central rapidities the T may still be 
antishadowed, but would be detected in the shadowing region with the muon spectrometer 



|24[. At the LHC, both J/ip and T production are shadowed. 

The shadowing and antishadowing features of Sa are shown explicitly in Fig. 2 where 
we compare S^ and Sj to the E772 data at 800 GeV on C, Ca, Fe, and W targets p3 . 
Recall that when Xf is increasing, x 2 is decreasing. The solid curves show S^, the dashed, 
S^ , and the short dashed, S'J. The difference between J/ip and ip' comes from their small 
mass difference. Antishadowing is clearly seen for the J/tp as x/ —>■ 0. Antishadowing at 
negative Xf is in contradiction with the data: clearly other effects are needed to account 
for this behavior. The effect of Q 2 evolution on S^ is apparent-the antishadowing peak is 



broader and closer to unity than the corresponding region of S 1 ^. The T enters the EMC 
region for Xf < 0, decreasing S^. While this trend agrees with the shape of the data, the 
magnitude is too large. 

Accounting for the primordial transverse momentum of the partons, qx, may reduce the 
magnitude of Sa since the transverse momentum spectrum may be broadened in a nucleus 
because of Fermi motion. If the cross sections are integrated over all px, the qx smearing 
will have no effect (the qx spectrum in a moving nucleon has the same normalization as for 
a stationary one). However, the data is taken over a finite range of px, suggesting that the 
broader qx smearing in the nucleus manifests itself in a reduced pr-integrated cross section 
relative to the nucleon, reducing Sa- We calculated the qx spectra in nuclei by smearing 
the primordial transverse momentum distribution in a free nucleon with nucleon momentum 
distributions from realistic many-body calculations |26|| . The primordial qx spectrum for 



a free nucleon was parametrized by a gaussian with (q T ) « lGeV 2 , the same form used 
to describe the low px Drell-Yan data in pp scattering f28fl . Similar values of (q T ) are 



also compatible with the low px data for J/ip production J29|. The convolution formula 
also involves longitudinal degrees of freedom and is therefore x dependent, producing an 
^/-dependent effect on Sa- A detailed description of our calculation will be shown in a 
more extended paper pTf . The nuclear qx smearing shown in Fig. 2(b) and in Fig. 2(d) 
for the J/ip (dotted curves) and in Fig. 2(d) for T (dot-dashed curve), reduces the effect 
of antishadowing and affects the shadowing part only marginally because of its increasing 
behavior with x (decreasing with xf). The agreement with data is substantially improved. 
In Fig. 3(a) we compare the model with the E772 Drell-Yan production data [[HJ in the 
W target as a check on the role of sea quark shadowing and the effect of Q 2 evolution. 
The model, including a small effect from qx smearing (dotted curve), indeed gives a good 
description of the data. We show a comparison with the NA3 J/ip production data from 



200 GeV proton interactions with H and Pt targets [29 in 3(b). The antishadowing region 



is broad with shadowing only for Xf > 0.4. The width of the antishadowing peak decreases 



with increasing energy since X2 passes from the shadowing region to the EMC regime over 
a smaller range of Xf. The antishadowing region is also shifted backward with Xf as the 
energy grows, indicated by a comparison of S 1 ^ in Fig. 2(d) and Fig. 3(b). Clearly nuclear 
shadowing alone, and even including the pt— integrated spectra, does not account for the 
observed suppression. Other effects become important and must be taken into account ||. 

Finally, we predict the shadowing of J/ip, T, and Drell-Yan production in Au+Au col- 
lisions at RHIC energy relative to pp production at the same energy in Fig. 3(c). Both the 
projectile and target distributions are shadowed. RHIC data will not be strongly affected 
by antishadowing, as expected from Fig. 1(c). Since only relatively small Xf values can be 
measured, features important at large Xf will not play a role. Thus shadowing will dominate 
the nuclear effects. We expect this behavior to be even more pronounced for LHC, where 
the range in x extends below x = 1CT 3 , Fig. 1(d). Here recombination effects should play an 
important role and could be tested accurately. Calculations including recombination effects 
at LHC energies are in progress and will be shown elsewhere [^7j . 

We have evaluated the effect of shadowing on heavy quarkonium production in current 
experiments and at RHIC. Because of their difference in mass, J/ip and T production follow 
different patterns. Antishadowing has been shown to be substantial at present energies, in 
particular, T production is within the antishadowing region. However the amount of shad- 
owing and antishadowing for T is smaller than for the J/ip because of the perturbative Q 2 
evolution of the shadowed distributions. Nuclear modifications of the primordial transverse 
momentum spectrum are an important addition in the low Pt range. 

Nuclear shadowing is an important feature of quarkonium production, yet at lower en- 
ergies it is not sufficient. Shadowing will play a major role in nucleus-nucleus collisions, 
suppressing both J/tp and T production significantly but according to different patterns in 
Xf, over an extrapolation from pp interactions. 
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Figure Captions 



Figure 1. A comparison of x 2 as a function of xj for J/ip (solid) and T (dashed) production, 
at (a) 200 GeV and (b) 800 GeV fixed target energies and (c) RHIC and (d) LHC ion collider 
energies. 

Figure 2. Our model calculations for S A (solid), S A (dashed), and S\ (dotted) compared to 
data from (a) C, (b) Ca, (c) Fe, and (d) W targets at 800 GeV ||19|| . The T data is shown on 
(d) only. The effect of q? broadening in the nucleus is shown in (b) and (d) for S A (dotted) 
and in (d) for T (dot-dashed). 



Figure 3. (a) A calculation of S A Y with the data from Ref. p3[ at 800 GeV on the W target 
(solid). The dotted curve shows the effect of qx broadening, (b) The J/ip production data 
of Ref. [^] on 2 H and Pt targets at 200 GeV compared with S A . (c) Predictions of S AA 
(solid), S AA (dashed), S AA (short-dashed), and S AA (dotted) for Au+Au collisions at RHIC. 
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